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Starch-water slurries are commonly used to study fracture dynamics. Drying starch-cakes benefit 
from being simple, economical, and reproducible systems, and have been used to model desiccation 
fracture in soils, thin film fracture in paint, and columnar joints in lava. In this paper, the physical 
properties of starch-water mixtures are studied, and used to interpret and develop a multiphase 
transport model of drying. Starch-cakes are observed to have a nonlinear elastic modulus, and 
a desiccation strain that is comparable to that generated by their maximum achievable capillary 
pressure. It is shown that a large material porosity is divided between pore spaces between starch 
grains, and pores within starch grains. This division of pore space leads to two distinct drying 
regimes, controlled by liquid and vapor transport of water, respectively. The relatively unique ability 
for drying starch to generate columnar fracture patterns is shown to be linked to the unusually strong 
separation of these two transport mechanisms. 

PACS numbers: 62.20.mm,47.56.-|-r,89.75.Kd 



I. INTRODUCTION 



grains themselves. 



Mixtures of starch grains and water have attracted at- 
tention recently as model materials for studying drying 
and fracture patterns, and in particular the formation of 
columnar joints [Ilia[3SlElEl[71lllll[ini[IIllia[ni 
[Mj [15] . When a starch- water slurry dries, it responds 
to desiccation in two distinct ways. Initially, the evap- 
orating slurry solidifies into a homogeneous starch-cake, 
which can be up to several cm thick, that dries uniformly 
throughout its thickness. As shown in Fig. [ija), stresses 
that occur in this phase drive the formation of a thin- 
film fracture pattern, similar to the craqulere patterns 
seen in pottery glazes or dried mud puddles [3l [TOl E] . 
Further evaporation leads to a form of directional drying, 
as a desiccation front is initiated at the drying surface, 
and intrudes into the starch-cake yij. This type of drying 
is associated with the development of columnar joints 
[T5] . such as those shown in Fig. [l|b). The resulting 
polygonal pattern is statistically identical to that found 
in columnar basalts, such as the Giant's Causeway [5], or 
Fingal's Cave, shown in Fig. [ijc). Interestingly, the col- 
umn size is also chosen by the same dimensionless scaling 
law in both cases |15j . It remains mysterious, however, 
why drying starch slurries should form columnar joints 
at all. Most other slurries or dispersions have only been 
seen to give rise to the more common craqulere patterns, 
when dried. 

Here, it is shown that the effective separation of the 
two very different drying mechanisms of liquid capillary 
flow and diffusive vapor transport is responsible for set- 
ting up the conditions of directional solidification that are 
necessary to produce columnar joints in a drying starch- 
cake. In the case of starch, this effective separation may 
arise from an unusually porous structure of the starch 
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FIG. 1: [Color online] Two types of cracks form in a drying 
starch-cake: (a) thin-film, or 'mud-crack' patterns, and (b) 
columnar joints. The thin-film crack spacing scales with film 
thickness [101 111], and is typically a few cm. The columnar 
joint crack spacing scales with the evaporation rate [121 114) . 
and is typically around a mm. (c) Columnar joints are com- 
mon in cooled lava, where they appear to be governed by 
similar scaling and organizational behavior to starch columns. 
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The stress state within a drying poroelastic solid is 
analogous to that of a cooling viscoelastic solid P^Bl ITT] . 
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This correspondence allows experimental results ob- 
tained from drying starch to be compared, quantitatively, 
with the thermal problem of columnar joint formation in 
lava [21 [TS] . Motivated by this well-known geophysical 
application, the dynamics of the columnar crack pattern 
have been studied experimentally in considerable detail 
[Illllll|71[HliaiIllIl[Tn]. starch suspensions can also 
act as a conceptual bridge between colloids and soils, as 
the size of starch grains is intermediate between those 
of most soils, and those of common colloids. There has 
been interest in studying water transport in starch as a 
soil analog [51 |B] , and to assist in explaining desiccation 
cracks in fine arid soils [IH]. The system has also at- 
tracted attention as a model test of simulations of crack 
formation [T31 [HI [501 (HI [12 , as the ordered columnar 
joints are expected to arise spontaneously from the un- 
derlying dynamics. In addition to quantiatively address- 
ing the physical origin of columnar joints in drying starch, 
it is hoped that the full description of drying and crack- 
ing of starch presented in this paper will encourage such 
physically based modeling, and foster the development of 
a common understanding of the formation of of thin-film 
fracture, desiccation cracks in soils, and columnar joints. 

The drying, and eventual fracture, of a slurry is a com- 
plex problem. A detailed understanding of the physical 
properties of the network of grains and pores is neces- 
sarily involved in any effort to model desiccation, and 
columnar jointing. Investigations were first made into 
several physical properties of desiccating starch slurries, 
in order to characterize how water may be transported 
either as a liquid or vapor, within an unsaturated starch- 
cake. The results presented here are then used to develop 
and test a dynamical model of water transport, and of 
the microscopic origins of stresses, within a desiccating 
starch-cake, initially composed of equal masses of corn- 
starch and water. Finally, this model is used to sug- 
gest why starches, unlike most other porous media, form 
columnar fracture patterns when they dry. 



II. PHYSICAL PROPERTIES OF 
STARCH- WATER MIXTURES 

For all experiments presented here, cornstarch (some- 
times known as cornflour) slurries were prepared by mix- 
ing equal masses of dry cornstarch (Canada Brand) and 
water in a glass container. Starch-cakes, up to 4 cm thick, 
were prepared by pouring the initial slurry into 60 mm 
radius flat-bottomed glass dishes, which were then placed 
under a pair of 250 W heat lamps, to dry. The materials 
and methods are identical to those reported in a number 
of studies of columnar jointing in starch [8[ [T2l ITS] . 

A. Particle size 

The properties of granular porous media are strongly 
affected by the size and geometry of their consitiuent par- 



ticles. Using a laser-diffraction particle sizer (Malvern 
Mastersizer S), the average particle radius of starch 
grains in water was found to be Rg — 8.2±0.1 /im, with 
the peak characterized by a full-width-half-maximum of 
3.2±0.1 /im. The starch grains do not swell when wetted 
- 8 tests were performed, from immediately after mixing 
dry starch with water, to 15 minutes later, with no vari- 
ation in size. 3 additional tests of grain size were done in 
air. Some clumping of grains in air was noticed, but the 
position of the main peak of the particle size distribution 
was not displaced from that measured in water. Particle 
size results were confirmed by the direct observation of 
a dilute sediment of starch in water, through a digital 
microscope. 



B. Particle surface character 

Electron microscopy was performed on a starch column 
selected from a fully dried columnar sample, and on dry 
starch grains dusted onto a sample holder. The samples 
were sputtered with a thin layer of carbon, and observed 
in a JEOL JSM-840 scanning electron microscope. When 
dried into a column, one surface of which is shown in 
Fig. |2ja), starch grains appear to be randomly close 
packed. On close inspection, as shown in Fig. [2](b), the 
grain surface is somewhat platy and broken up, and the 
grains are seen to contain many pits. When wetted these 
pits can potentially fill with water, which will need to 
be considered when estimating the porosity of any grain 
packing. The large surface area of the pits, and their 
small diameter, may also account for some of the strong 
hydrophilic qualities of corn starch. 

If a starch suspension is allowed to dry and form colum- 
nar joints, the starch grains are damaged on the micro- 
scopic level. If a sample is re-wetted, mixed, and re- 
dried, the resulting columns are usually larger and more 
irregular than those resulting from the first drying. This 
damage is noticeable at grain-grain contacts in desiccated 
samples. As shown in Fig. 2(c), these contacts involved 
the compression of one grain into the other, or the flatten- 
ing out of sections of both grains. It is likely that these 
'kissing' contacts were created by the action of water- 
bridges joining the neighboring grains. Such bridges cre- 
ate a strong capillary force, pulling adjacent grains to- 
gether. Any strength of the final, desiccated columnar 
state is presumably due to the surface (van der Waals) 
forces acting across these contacts, as it is in powder 
compacts [23]. 

C. Riieology 

Granular solids are known to have highly non-linear 
rheologies. Since stress is only transmitted via particle- 
particle contacts, which behave as Hertzian springs |24j . 
a non-linear elastic model is appropriate. The simplest 
such model to approximate a granular solid is that of 3rd 
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FIG. 2: Scanning electron microscope images of starch, (s) 
Starch observed on the surface of a columnar joint are ran- 
domly, and closely packed, (b) The starch grains are roughly 
spherical, and display both a platy, broken surface, and a 
highly porous sub-structure, (c) Individual grain-grain con- 
tacts show deformation of the grains. 

order elasticity. This model was originally developed by 
Landau to describe anharmonic elastic effects in crystals 
[25j . but has been successfully applied to granular media 
in a wide variety of engineering and geophysical contexts 
[55] . In its general form this model adds the next-lowest 
(cubic) order terms to the elastic energy density. In the 
absence of shear, the stress tensor 

a,, = Eeu + Gel, (1) 

where E is the Young's modulus, and C is one of three 
general Landau moduli, which describes the rate at which 
the stiffness of a body increases with increasing strain 
|25j . The effective Young's modulus E' = an /en was 
measured in a 2 mm thick, partially dried starch-cake by 
compressing and releasing the sample under a 40 mm di- 
ameter, parallel plate rheometer headpiece (TA Instru- 
ments RAIOOO), while observing the sample thickness 
and normal force response. Experiments were performed 
over a range of compression and relaxation speeds, from 
0.1-4 /im/s, and there was no observable effect of the 
strain rate on the normal force response. As shown in 
Fig. [sjja), E, the elastic modulus at zero strain, was 
found to be consistent with zero, and the Landau mod- 
ulus C=16±l MPa. Tests at moisture concentrations of 
0.3 and 0.4 g/cm'^ showed no difference in results. 

An unusual feature of some granular solids is their lack 
of a well-defined stress relaxation time. They can display 
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FIG. 3; Drying starch exhibits highly non-lilnear rheology. 
(a) The effective Young's modulus in a 2 mm thick starch-cake 
of C = 0.3 g/cm^ is linearly dependent on the vertical strain 
(tests at strain rates of 0.1 /im/s (black), 0.5 /im/s (grey), 
and 2 /xm/s (white) are shown), (b) In a starch-cake the 
logarithmic relaxation of compressional stress (black circles, 
left axis, two independent trials shown) and shear-strain (grey 
circles, right axis) continues for a long time after a stress or 
shear is applied. 



a slow dynamical relaxation of an applied stress, that 
is logarithmic in time [26]. This behavior belongs to a 
class of relaxation mechanisms shared by a number of su- 
perficially unconnected systems in energetically jammed 
states, such as of the magnetization of spin glasses, or 
the dc susceptibility of granular magnetic media [27] . 
The causes of these slow dynamics are not particularly 
well understood, although Pauchard et al. have recently 
shown how the slow viscous deformation of nano-particles 
can relax stress [28]. The stress-relaxation of starch-cakes 
was tested by compressing a cake with a fixed strain, and 
monitoring the responding normal force over time. Al- 
ternatively, a fixed shear stress was applied for 5 minutes, 
and the relaxation of the shear strain was monitored, af- 
ter the stress was released. As shown in Fig. |3jb), all 
starch-cakes showed a stress or residual shear that de- 
cayed logarithmically in time. These results imply that 
stress relaxation will only play a small role of the rheo- 
logical response of a starch-cake to stress; over the time- 
scale of a desiccation experiment, lasting of order 10^-10^ 
s, these relaxation mechanisms could only reduce an im- 
posed stress by approximately 15-20%. Since the starch- 
cake remains irreversibly deformed after drying, however, 
some form of permanent particle deformation must occur. 
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three-phase (starch- water-air) mixture, is described by 

{1 - Ps/pw)Ps + C, C>cj)pw, 



P = 



(2) 



where ps is the bulk density of dry starch, p^, is the den- 
sity of water, and (j) is the porosity of the packed starch- 
cake. A fit of Eqn. |2]to the density of starch mixtures, 
as shown in Fig. 4][a), yields ps = 1.57±0.01 g/cm'^, and 
(j) = 0.56±0.01, and is in reasonable agreement with a 
previous measurement of the porosity of starch of (j) ~ 
0.51 [6 . 

Below a moisture concentration of approximately 0.3 
g/cm^, an unstrained three-phase density model is no 
longer valid, as the starch-cake responds to further drying 
by shrinking. The average vertical strain in these sam- 
ples was estimated by comparing the actual thickness of 
the starch-cakes with the expected unstrained thickness. 
As shown in Fig. Qb), the onset of strain is well cor- 
related with the appearance of columnar jointing in the 
starch-cake. The total strain measured is quite large, ap- 
proximately 0.1. Although some deformation must occur 
before the starch-cake dries to C = 0.3 g/cm"^, as primary 
cracks open during this phase of drying, the methods pre- 
sented here could not detect the small strains involved. 



FIG. 4: Density changes and strain develops during the dry- 
ing of a starch slurry, (a) the measured density dependence 
(open circles) of starch-water mixtures on C is well mod- 
eled by a non-interacting mixture of water, air, and starch 
(dashed curve), (b) below C = 0.3 g/cm"^ the density devi- 
ates from this curve due to desiccation-induced strain. The 
grayed region indicates samples in which columnar joints were 
observed, but where joints did not reach the base of the starch- 
cake. 



This was observed to be the case at the kissing contacts 
shown in Fig. |2jc). 

D. Density 

The density of drying starch slurries was studied, in or- 
der to measure the porosity of a drying starch-cake, and 
to observe the onset and development of strain in the dry- 
ing mixture. Liquid starch-water mixtures were added 
to a 100 ml volumetric flask and weighed to determine 
their density. An oscillating U-tube densitometer (Anton 
Paar DMA 500) was also used to measure fluid density, 
although for dilute samples the settling of starch out of 
suspension affected these measurements. The densities of 
solid mixtures were calculated from measurements of the 
thickness, diameter, and mass of partially dried starch- 
cakes. 

If no tensile strain occurs during drying, then the de- 
pendence of the density p of a starch- water mixture, on 
its water concentration C, is well flt by a simple model of 
non-interacting phases. This mixture density, p, across a 
transition from a two-phase (starch-water) mixture to a 



E. Permeability 

The transport of liquid water in a drying starch-cake 
occurs by flow through a porous medium. In these con- 
ditions, Darcy's law 



(3) 



describes the mass flow rate per unit area, q; , of a fluid of 
viscosity p and density p that will pass through a porous 
body of intrinsic permeability n, under a pressure gra- 
dient VP. The permeability was measured in order to 
characterize the efficacy of liquid transport during dry- 
ing. It is also used as a probe of the inter-particle poros- 
ity, (f)p, which may differ from the total porosity = 0.56, 
due to the porous structure of the starch grains. 

A value of k =(4.5±1.5)x10^^'* m^ was measured in 
saturated starch samples, at 24°C. A filter was glued to 
the bottom end of an open glass tube, and starch slurry 
was poured in from the top and allowed to settle. Wa- 
ter was then added above the starch plug, and the flow 
rate through the starch was measured. There was no 
dependence of the permeability on sample thickness or 
hydraulic head, and measurements were stable over long 
periods of time; one test was run for 80 hours, with no 
significant drift in the measured k. 

The Carmen-Kozeny equation 



45(1 - 



"12 



(4) 



is a commonly used semi-empirical model of the perme- 
ability of a porous medium ^29. ,30) . Using a particle 



5 



100 



10 



0.1 





1 ' 1 


' 1 ' 1 - 


M 




^ K>\ 






























m 












m ^ 

m 




m 

m 


« K 






m 

w 

1,1 








0.1 



0.01 



0.05 0.1 0.15 0.2 0.25 
Water concentration (g/cm^) 



0.3 



FIG. 5: The water potential (open circles) in partially dried 
starch-cakes increases as the cake dries. The radius of curva- 
ture (filled circles) of water pockets in equilibrium with this 
potential drops dramatically as the cake dries. Values above 
C =0.25 g/cm^ are not within detectable limits. 



radius of 8.2 /xm, and the measured permeability, Eqn. 
[4] predicts that (ftp = 0.26±0.03, which suggests that the 
inter-particle porosity is much lower than the total poros- 
ity. 



representative of the interparticle spaces, which are ap- 
proximately an order of magnitude smaller than Rg [32] . 
Below C = 0.25 g/cm"^ the water potential rose rapidly 
as the starch-cakes dried, indicating that any remaining 
water is trapped in a distribution of much smaller pores. 



G. Fluid phase 

25 g of starch were mixed with 500 ml of water for 
1 hour using a stir bar. The resulting suspension was 
left to settle overnight, and decanted. Measurements 
of the density of the supernate were indistinguishable 
from measurements of the density of tap-water. Simi- 
larly, measurements of viscosity agreed with that of tap- 
water, below the gelation temperature Tg = 62° C. Above 
Tg the supernate's viscosity increased to approximately 
twice that of water, indicating a trace of dissolved amy- 
lose in solution ^33j. These results confirm that, as long 
as experiments are kept below Tg, the fluid phase of a 
starch-cake can be treated as pure water. 



III. DRYING DYNAMICS 



F. Water potential 

In addition to fluid flow, water can be transported 
within a drying starch-cake by vapor diffusion. These 
two mechanisms are in fact linked, as water can change 
phase. The equilibrium vapor pressure of water, Py, over 
a liquid surface varies according to the Kelvin equation 

P. - P.a.e^^"/«^, (5) 

where R is the gas constant, Vm is the molar volume of 
water, P^at is the equilibrium vapor pressure of water at 
temperature T over a flat liquid surface , and the pressure 
drop 



is caused by a surface tension F acting across an air- 
liquid interface with a radius of curvature R^, |31| . For 
the concave menisci in a porous medium, Py is slightly 
reduced from Psat- 

The water retention curve (also known as the soil wa- 
ter characteristic curve) describes how P changes with 
C, and can also be used to estimate the pore-size dis- 
tribution of a material. Shown in Fig. |5] points along 
this curve were measured using a Decagon WPT4 water 
potentiometer. 1-2 mm thick starch-cakes were dried to 
a desired C, placed in the test cell, and allowed to equi- 
librate. P was then measured by observing the relative 
humidity in the test cell. Signiflcantly above C = 0.25 
g/cm'^ the water potential was consistently below a de- 
tection limit of ~0.5 MPa. This is expected when R^ is 



The drying of a starch slurry occurs in three stages, as 
has been widely recognized [3 [HI H [TUl HU El US] . These 
are demonstrated in Fig. |6] Each phase is contoUed by 
a different water transport process. 

Since most desiccation experiments begin with an ex- 
cess of the fluid phase, so that the mixture can be ma- 
nipulated as a liquid, the water content is initially higher 
than about 0.6 g/cm.^^. The first stage of drying occurs as 
the particles settle in a water column, while water evap- 
orates from the surface of a pool of supernatant liquid. 
The surface evaporation rate is constant, and has been 
used to characterize the relative drying rate of starch- 
cakes exposed to different drying conditions [3 [T^] . 

In order to describe the dynamics of the later stages 
of drying, the water concentration C{z,t) was sampled 
at different vertical heights z, measured from the base of 
the starch-cake, and times t in a series of 17 identically 
prepared samples, each of which dried according to the 
evaporation rate shown in Fig. [6j Replicates were re- 
moved from the drying apparatus after different times t 
and sliced into thin layers that were then weighed, fully 
dried, and reweighed in order to measure the water con- 
tent. The position of the columnar fracture front in each 
sample was also recorded. The resulting moisture con- 
centration field, shown in Fig. |7] sampled C at roughly 
3 hour intervals over 1-2 mm thick layers. 

When a starch-cake dries to C = (j)pw, air begins to 
infiltrate into the voids between starch grains, initiat- 
ing the so-called funicular regime of drying [53]. This 
regime, also sometimes called the 'first falling rate' pe- 
riod, is dominated by fiuid fiow through capillary bridges 
[33] . There is little to no measurable strain during this 
regime, although primary cracks do appear. By this time. 
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FIG. 6: The drying of a starch suspension, initially mixed 
with equal masses of water and starch, can be divided into 
three stages. From the left to the right: (1) starch grains 
sediment in a water column, with surface evaporation; (2) a 
rigid mixture of starch, water, and air maintains a high surface 
evaporation, as capillary flow transports water to the drying 
surface; (3) once the capillary network becomes disconnected, 
the slower process of vapor transport completes drying. 
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FIG. 7; [Color online] Moisture concentration C{z,t) mea- 
sured in a series of identically prepared, partially dried starch- 
cakes. Experiments were conducted to sample C approxi- 
mately every 3 h, and every 2 mm - the displayed field is 
linearly interpolated between data, and contours are drawn 
every 0.05 g/cm^. The fracture front position (open circles) 
remains at C/ = 0.20±0.01 g/cm^, slightly trailing the begin- 
ning of the drying front. 



the starch-cake is solid, has the consistency of moist chalk 
or clay, and is very weak and brittle. The cake, however, 
still contains water equivalent to up to 56% of its vol- 
ume. A random close packing of spheres, which might be 
expected to roughly characterize the grain structure, has 
a porosity of only 36%. In order to realize such a large 
porosity, it appears likely that the water in a starch-cake 
is distributed between two distinct but connected reser- 
voirs: the gaps left between the packed grains {4>p), and 
the smaller spaces within individual grains {4'g)- 

As evaporation progresses, larger pores will empty be- 



fore smaller ones in equilibrium with them. Thus, if there 
is a broad distribution of pore types, once most of the wa- 
ter has been removed from the pore space between grains, 
the remainder will be confined to small water bridges be- 
tween particles, and any smaller pores within particles. 

As its water content falls past 0.30 g/cm'^, there is 
a sudden change in the drying behavior of the starch- 
cake. In similar cases, such a drop in the drying rate 
indicates a transition from a network of connected capil- 
lary bridges, to an unconnected network [131 [34]. At the 
onset of this 'pendular' regime, capillary flow abruptly 
ceases, and moisture can be transported only through 
the much slower means of vapor diffusion. For a random 
close-packing of spheres, the funicular-pendular transi- 
tion should occur at around C = 0.05 g/cm^ [3S]. In 
order to explain the high value of C observed at this 
transition, 0^ must contribute an additional porosity of 
approximately 0.25. 

During the pendular phase, a drying front forms, and 
passes through the sample, from top to bottom. At 
this front, as shown in Fig. [6] there is a clearly delimi- 
nated boundary between nearly uniformly moist unfrac- 
tured starch, and starch which has given way to colum- 
nar joints. In contrast to a typical diffusive front, the 
curvature d^C/dz^ is positive above the fracture front 
|12j . Similar sharp drying fronts have been seen in the 
directional drying of 2D suspensions ^29| , and have been 
imaged in drying starch slurries using nuclear magnetic 
resonance techniques [9 . Such a front is expected when 
there is a shift from flow-limited to diffusion-limited wa- 
ter transport [1^01^ . 



IV. WATER TRANSPORT MODEL 

In order to quantify the remarks made in the previous 
section, a model of water transport is briefly outlined 
here. The basic model was recently proposed by Nishi- 
moto et al. [13 , and further details are given there. The 
model presented here has been modifled to account for 
the implications of pore space within starch grains. 

Fluid flow in a porous structure is governed by Darcy's 
law, Eqn. [3] with a permeability k(C) that depends on 
the degree of saturation. Expressed in terms of a concen- 
tration gradient, this implies that the fluid flux 



k{c)p^ dP 

dC 



(7) 



Here, variations in the driving pressure P arise from the 
dependance of the radius of curvature of the capillary 
bridges connecting particles, on C . 

The transport of water vapor is diffusive. Since the 
equilibrium vapor concentration is proportional to the 
vapor pressure P„, the Kelvin equation, Eqn. |5j can be 
used to show that 



q« = 



satVn 

tRT 



(8) 
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FIG. 8: The variation of the diflusivity of water within a 
drying starch cake can be measured by either considering the 
evolution of the conce ntra tion (Eqn. 

water potential (Eqn. 10 filled circles) . The water potential 
: 0.3 cmVs and Csat = 5.1x 10"^ g/cm^* 
40° C of the starch- 
= 0.56, and that K, 
dC/dz is small, and 



competing empirical models, but demonstrate how the 
difFusivity increases rapidly between C = 0.25 and 0.3 
g/cm"^. This suggests that at this water concentration, 
the starch-cake is indeed reaching the pendular-funicular 
transition, with an estimated (pg — 0.25. 

In the pendular regime, when the fluid network is un- 
connected, it is assumed that all moisture transport oc- 
curs within the vapor phase, and that k. — Q. Under 
this assumption, as shown in Fig. |8] the water poten- 
tial P(C) can be used to estimate the diflusivity D{C) 
directly, using Eqn. |10| This interpretation of the wa- 
ter potential is compared to the solution of Eqn. [TT] in 
Fig. |8] and accurately captures both the beginning of 
an rapid rise in D(C) at low moisture concentrations, 
and a plateau in D(C) at intermediate concentrations. 
The small discrepancy, approximately a factor of two, 

may 



open circles), or the between D{C) as calculated from Eqns. ITo] and 
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data assume that Do ■ 
are taken at a typical temperature T 
cake, that r = \/2, = 1 g/cm^ 
= 18 cm^/mol. Above C = 0.3 g/cm^, 



Eqn. [TT] cannot be accurately calculated. 



where Csat is the saturated concentration of water in air. 
Do is the free diflusion constant of water vapor in air, 
and r is the tortuosity of the network through which the 
gas must diffuse [T^ . 

By combining the liquid and vapor fluxes with a mass 
conservation equation, one obtains a nonlinear diflusion 
equation 



dC 

1h 



= -V-(q, + q„) = V-p(C)VC), 



(9) 



which has previously been used to describe mass trans- 
port in a drying slurry [9l [121 1131 US]. In this case, the 
diflusivity 



D = 



DqCsatVrn ,^ C p^^ / Hj. ^ k{C) Pyj \dP 

tRT M / dC 



Alternatively, by integrating Eqn. [9] and by assuming 
that the concentration varies only over the vertical posi- 
tion < z < /i, D{C) can be calculated from the mea- 
surements of the concentration field presented in Fig. [7] 

mm, 



D{C{z,t)) = 



dC 

In 



dz' 



dC 
dz 



(11) 



Evaluating Eqn. [9] requires considerable detailed 
knowledge of the system under study. In the funicular 
regime, there exist a number of semi-empirical models of 
capillary transport [T31 [3U [33], which typically assume 
either a power law |13j or exponential |36| increase in the 
diflusivity with increasing water concentration. In Fig. [8] 
the measured concentration field has been used to calcu- 
late D{C). These data do not discriminate between the 



be attributable to diflusive transport through the starch 
grains, to local transport of water within pores, or to 
temperature gradients within the starch-cake (Csat is 
strongly temperature-dependant), which are not included 
in the model. The reasonable agreement between the 
two methods implies that Eqn. [8] accurately describes 
a vapor-controlled transport of moisture in the drying 
starch-cake, below C ~ 0.20 g/cm^. 

The transition between liquid and vapor transport 
mechanisms is characterized by a minimum in D(C), 
which creates a bottleneck in the water transport to 
the drying surface. The existence of some minimum in 
D{C) has been described [T2[T31[IS], but a broad, deep 
minimum appears to be necessary for creating the sharp 
drying front associated with the columnar fracture front 
|15| . The depth and position of this minimum depends 
on the water storage capacity of the starch grains. In 
the absence of appreciable pore space within grains, the 
pendular regime would be confined to C < 0.05 g/cm^ 
[55] . Under these conditions, Nishimoto et al. used a 
semi-empirical model, adapted from soil mechanics, to 
evaluate Eqn. |11[ They predicted a minimum in the dif- 
fusivity that was one-to-two orders of magnitude larger 
than observed [13J. In this case, the resulting moisture 
distribution, C(z, t), would be considerably more difluse. 



V. FRACTURE MECHANICS 

There are two types of fracture that have been studied 
in drying starch-cakes. While the cake is in the pendular 
regime there is eflectively no variation of water concen- 
tration on position, and first-generation (also known as 
primary) cracks appear that penetrate the full depth of 
the drying cake. Later, during the funicular regime, dry- 
ing is essentially confined to a thin interface between the 
wet and dry regions of the starch-cake. Columnar joints 
form during this phase, with the crack tips following the 
drying front. 

The plumose structure on the primary cracks show that 
they initiate at the upper drying surface, and quickly 
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propagate vertically through the entire sample [5]. The 
resulting disordered fracture networks scale like those of 
ordinary thin films, and have been studied as examples 
thereof [TOJ [TT]. In particular, the average spacing be- 
tween first generation cracks increases roughly linearly 
with sample depth, and these cracks typically meet at 
90° junctions [11 [TT]. 

As it enters the funicular regime, a starch-cake begins 
to dry from the top down. Water is confined in uncon- 
nected reservoirs within particles, and in inter-particle 
capillary bridges. During the final phase of drying, a to- 
tal strain of up to 0.1 develops, which, according to Eqn. 
[l] is equivalent to putting the sample under a uniform bi- 
axial stress of up to 160 kPa (the vertical stress can be 
relieved by lowering the drying surface). Internally, this 
strain is accommodated by the deformation of grains into 
each other, as shown in Fig. |2] through the process of 
liquid sintering [37]. This measured stress is in reason- 
able agreement with the expected maximum attainable 
capillary pressure, typically taken to be of order I2"f/Rg 
= 110 kPa [Ml [31]. 

A fine-scale fracture pattern develops at the drying 
surface at the start of the funicular regime, and further 
drying propagates the crack tips into the sample. This 
leads to a moving fracture front that tracks a clearly de- 
lineated boundary between relatively soft, moist starch, 
and harder dry starch, which has given way to columnar 
jointing. These two layers come apart easily when a par- 
tially dried sample is removed from its dish, as all the 
fracture tips are confined to a thin interface. 

The directional drying of an unsaturated starch-cake 
is related to the directional drying of thin colloidal lay- 
ers [23 ISSl HO]- In these experiments, colloidal disper- 
sions are typically confined within a small cell, and al- 
lowed to dry from one exposed face. Stress builds up in 
the drying body, in response to a gradient in pore pres- 
sure. Eventually, an array of cracks form at the exposed 
edge, and slowly advance into the solidifying network. 
Like columnar joints, these cracks are evenly spaced out, 
and advance perpendicular to a drying front. However, 
unlike joints in drying starch, it should be pointed out 
that cracks in colloidal dispersions occur while the drying 
body is still saturated with water [25] . 

Throughout unsaturated drying, the surface tension of 
the water bridges between grains gives a partially dried 
starch-cake much of its cohesiveness, and acts to pull 
the starch grains together. The bridge-grain structure 
has been modeled as a spring whose equilibrium posi- 
tion varies with C, and which breaks when stretched be- 
yond finite displacement [13i il9. ,38, . The compaction of 
grains into each other leads to the development of cracks, 
as any set of particle-particle contacts will contain some 
bonds which are weaker than others. By allowing such 
bonds to yield, the surrounding network of bridges re- 
laxes. These assumptions give rise to fracture patterns 
that agree qualitatively with experiments, although there 
has been some difficulty in studying model systems suffi- 
ciently large enough to be checked quantitatively against 



experimental observation [T31 [T5J |3S]. Future develop- 
ments of this type of discrete-element model would ide- 
ally couple the water transport dynamics described above 
with the evolution equations for the springs. 

Columnar joints have been reported in many thermoe- 
lastic systems, such as lava, smelter slag, and optical 
glass [41 J. However, starch remains the only known sys- 
tem in which columnar joints have been studied in de- 
tail. Indeed, this unusual feature of dried starch-cakes 
has been independently rediscovered on several occasions 
[TJ 135] [55] . It has recently been shown that the diameter 
of columnar joints scales with the thickness of the drying 
front |15j . which is in turn controlled by the cross-over 
in drying dynamics from flow-dominated transport, to 
vapor-dominated transport. In drying starches, this front 
is typically a few mm wide ^15j , and columnar joints are 
a few mm in diameter. 

It has been argued that columnar jointing occurs when 
the length scale set by the minimum diffusivity, and a 
characteristic drying time, is smaller than the thickness 
of the drying layer |15j . In the previous section, it was 
demonstrated that the sudden transition between wet 
and dry starch arises from a broad, deep minimum in the 
diffusivity of moisture throughout the starchcake, which 
has been interpreted to arise from a large internal poros- 
ity, within starch grains. In the absence of these pores, 
the minimum diffusivity has been estimated to be sig- 
nificantly larger [T3], and the drying front would have 
a characteristic width of many cm. This suggests why, 
in most natural or laboratory situations, given that most 
soils do not share the observed porous structure of starch 
grains, the evolution of desiccation fracture patterns typ- 
ically halts at the mud-crack stage. In other cases it 
may be possible to produce columnar joints, if the dry- 
ing body is large enough. The experiments of Gauthier 
et al. |40l , who observed the slow propagation of two per- 
pendicular cracks, meeting at a vertex, in colloidal silica 
confined within a 1 mm diameter tube, suggest that 3- 
dimensional crack patterns such as columnar joints may 
be experimentally accessible in other drying systems of 
modest scale, for example. 



VI. CONCLUSION 

Starch-water mixtures shrink and crack when dried. 
The simplicity of this system, and its use as a model for 
studying columnar joints, have encouraged considerable 
interest in the drying and cracking behavior of starch- 
cakes [I][l[31ll[51[Sl|7]|H]ini|Tn]III]IIllI3[IllISi. 

The results presented here have shown, quantitatively, 
how directional drying occurs in starch-water mixtures. 
An effective separation of vapor and liquid water trans- 
port creates a bottleneck in the drying behavior, leading 
to this type of drying. 

In order to evaluate how water is transported through 
a starch-cake, several properties of starch were measured. 
An unusually large cake porosity was observed, which was 
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argued to be divided roughly equally between pore space 
within grains, and pore space between grains. By adapt- 
ing a model of moisture transport, it was shown how this 
division affects the transition between flow-dominated, 
and vapor-dominated water transport, causing the liquid 
transport to become ineffective in circumstances where 
considerable water remains in the body. A broad, deep 
minimum in the effective diffusivity of moisture followed, 
and lead to a sharp, propagating drying front. The total 
strain associated with this front is comparable to that 
predicted by capillary effects between grains. Such a 
sharp shrinkage front appears to be necessary to generate 
columnar joints |15] . The relatively unique behavior of 
drying starch slurries to give rise to columnar joints can 



thus be seen to arise directly from the unusual porous 
structure of starch grains, and the mechanics of moisture 
transport. 
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